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The motion of a stone skimming over a water surface is considered. A simplified description of 
the collisional process of the stone with water is proposed. The maximum number of bounces is 
estimated by considering both the slowing down of the stone and its angular stability. The conditions 
for a successful throw are discussed. 



I. INTRODUCTION 

Nearly everyone has tried to throw a stone on a lake 
and count the number of bounces the stone was able to 
make. Of course the more, the better.El Our intuition 
gives us some empirical rules for the best throw: the best 
stones are flat and rather circular; one has to throw them 
rather fast and with a small angle with the water surface; 
a small kick is given with a finger to give the stone a spin. 
Of course these rules can be understood using the laws of 
physics: the crucial part of the motion is the collisional 
process of the stone with the water surface. The water 
surface exerts a reaction (lift) force on the stone, allowing 
it to rebound. This process is quite complex because it 
involvgs-the description of the flow around the immersed 
stone.tlt!l Some energy is also dissipated during a collision, 
so that after a few rebounds, the initial kinetic energy of 
the stone is fully dissipated and the stone sinks. 

The purpose of this paper is to propose a simplified 
description of the bouncing process of a stone on water, 
in order to estimate the maximum number of bounces 
performed by the stone. This problem provides an enter- 
taining exercise for undergraduate students, with simple 
explanations for empirical laws that almost everyone has 
experienced. 



II. BASIC ASSUMPTIONS 

Consider a flat stone, with a small thickness and a mass 
M. The stone is thrown over a flat water surface. The 
angle between the stone surface and the water plane is 
8. A schematic view of the collisional process is shown 
in Fig. [|. The velocity V is assumed to lie in a sym- 
metry plane of the stone (the plane of the paper). The 
difficult part of the problem is, of course, to model the 
reaction force due to the water, which results from the 
flow around the stone during the stone- water contact. It 
is not the aim of this paper to give a detailed descrip- 
tion of the fluid flow around the colliding stone. Rather 
I shall use a simplified description of the force keeping 
only the main ingredients of the problem. First, the ve- 
locity V of the stone is expected to be (at least initially) 
the order of a few meters per second. For a stone with 
a characteristic size a of the order of a few centimeters, 
the Reynolds number, defined as Tie = Va/v, with v the 
kinematic viscosity (y ~ 10~ 6 m 2 s _1 for water), is of or- 



der Tie ~ 10 5 , that is, much larger than unity.Q In this 
(inertial) regime, the force due to the water on the stone 
is expected on dimensional grounds to be quadratic in 
the velocity and proportional to the apparent surface of 
the moving object and the mass density of the fluid.El Be- 
cause the stone is only partially immersed in water during 
the collisional process, we expect the force to be propor- 
tional to the immersed surface (see Fig. [I]). The force 
can be adequately decomposed into a component along 
the direction of the stone (that is, along t, see Fig. fi]) 
and a component perpendicular to it (that is, along n). 
The latter corresponds to the lift component of the force, 
and the former corresponds to a friction component (of 
water along the object). I write the reaction force due to 
water, F, as: 

F = ^C lPw V 2 S im n + ^C fPw V 2 S im t , (1) 

where C\ and Cf are the lift and friction coefficients, p w is 
the mass density of water, Si m is the area of the immersed 
surface, and n is the unit vector normal to the stone (see 
Fig. [j]). Note that in general, both Cj and Cf are func- 
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FIG. 1: Schematic view of the collisional process of a flat 
stone encountering a water surface. The stone has a velocity 
V, with an incidence angle (3, while 6 is the tilt angle of the 
stone. The immersed area Sim represents the area of the stone 
in contact with the water surface. The depth of the immersed 
edge is 2. 

tions of the tilt angle 9 and incidence angle (3, defined as 
the angle between velocity V and the horizontal. In the 
simplified analysis I will assume that both Ci and Cf am 
constant and independent of tilt and incidence angles.H 
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This assumption is not a strong one because ricochets 
are generally performed with a small tilt angle, 9, and a 
small incidence angle, (3. If one denotes the initial com- 
ponents of the incident velocity by V x q and V z q (parallel 
and perpendicular to the water surface, respectively), the 
latter assumption amounts to V z q <C V x q- 

We expect the lift force to be maximum when the ob- 
ject is only partially immersed due to the lack of symme- 
try between the two sides of the stone. Therefore, if the 
object reaches a depth such that it becomes completely 
immersed, the lift force would be greatly diminished and 
would probably not be able to sustain the weight of the 
stone anymore. For simplicity, I will assume that the 
lift force vanishes for completely immersed objects. The 
model for the force in Eq. (|l|) is crude, but it is expected 
to capture the main physical ingredients of the stone- 
water interaction. It might fail for lower stone velocities 
or larger incidence angles, where a bulge of water could 
be created and affect the lift and friction forces on the 
stone.El However, in this case it is expected that the stone 
will be strongly destabilized during the collision process 
and perform only a very small number of bounces. We 
will restrict ourselves to large initial velocities and small 
incidence angles, such that the number of bounces is suf- 
ficiently large. 



and it might fail in the last few rebounds of a stone skip 
sequence. 

With this approximation, Eq. (2b) decouples from 
Eq. (|2a|) . I thus first focus the discussion on the equation 
for the height z, which is the height of the immersed edge 
(see Fig. p. Note that the equation for z is equivalent 
to the equation of the center of mass position, Eq. ( pb| ) 
because 9 is assumed to be constant (see Sec. |v| for a de- 
tailed discussion of this point). Hence, we may identify 
V z with dz/dt and Eq. ( j2b| ) yields a closed equation for 
the height z. 



IV. COLLISIONAL PROCESS 

To solve Eq. (|l]) we need to prescribe the z depen- 
dence of the immersed area 5j m . This quantity depends 
on the precise shape of the stone. A n atural choice is 
circular, which I will treat in Sec. IVB. However, it is 



enlighting to first consider a square shape; this shape 
greatly simplifies the mathematics and already contains 
the basic mechanisms involved. 



A. A Square Stone 



III. EQUATIONS OF MOTION 

Consider the collisional process, that is, the time dur- 
ing which the stone is partially immersed in water. I 
will assume in this section that the incidence angle 9 be- 
tween the stone and the water surface is constant during 
the collisional process. The validity of this assumption is 
considered in detail in Sec. fv|. The origin of time, t = 0, 
corresponds to the instant when the edge of the stone 
reaches the water surface. During the collisional process, 
the equations of motion for the center of mass velocity 
are 

M ^T = -lpwV 2 S im {Cism9 + C f cos9) (2a) 
= -Mg + ^p w V 2 S im (Cicos9-C fS m9X?b) 

with V 2 — V 2 + V 2 and g is the acceleration due to 
gravity. Note that in Eq. (^) the area Si m depends on 
the immersed depth and thus varies during the collisional 
process. 

Equation (|J) is non-linear due to the V 2 terms on the 
right-hand side, but also due to the dependence of the 
immersed area, Si m , on the height z. However, we can 
propose a simple approximation scheme: the magnitude 
of the velocity, V, is not expected to be strongly affected 
by the collision process (as I shall show in Sec. VI). I 
thus make the approximation that V 2 ~ V 2 Q + V 2 Q ~ V 2 Q 
on the right-hand side of Eq. (^) . The validity of this as- 
sumption requires a sufficiently high initial velocity, V x q, 



In this case, the immersed area is simply S- lm = 
a|z|/sin# (see Fig. Q), with a the length of one edge of 
the stone. The equation for z thus becomes 



= -Mg ~ i Pw V 2 C— . (8i 



az 
sin( 



where C = Ci cos 9 — Cf sin 9 ~ Cj, and I have used 
\z\ = —z (z < 0). We define the characteristic frequency 
wo as 



2 _ Cp w V 2 a 



and rewrite Eq. (||) as 

fz 
dt 2 



2Msin( 



+ uj z = -g . 



(4) 



(5) 



With the initial conditions at t = (first contact with 
water), z — and i = V z o < 0, the solution of Eq. (||) is 

u\ 9.9 , , Vzo ■ , ,„\ 
z(t) = 2 "I — 2 COSUJ ot H smwoi . (d) 

Equation (|^) characterizes the collisional process of the 
stone with water. After a collision time t co u defined 
by the condition z(t co \i) — (t co u ~ the stone 

emerges totally from the water surface. It is easy to show 
that the maximal depth attained by the stone during the 
collision is 
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As discussed in Sec. |, the stone will rebound if it stays 
only partially immersed during the collision. The re- 
bound condition can be written as |z max | < asin#. If 
we use Eqs. (0) and (jj), this condition can be written 
after some straightforward calculations as 



V x0 >V C = 



iMg 



1 _ _2 tan 2 f3M 
a 3 Cp w sin 9 



(8) 



where the incidence angle j3 is defined as V z q/V x q — 
tan (3. Therefore, we obtain a minimum critical veloc- 
ity for skimming. Using the typical values, M = 0.1kg, 
L = 0.1m, Ci w C f w 1, p w = 1000 kgm~ 3 , and 
P ~ ~ 10°, we obtain V c ~ 0.71m s" 1 - lms" 1 . 

The physical meaning of this condition is clear: it sim- 
ply expresses the fact that the lift force \C p w V 2 a 2 has 
to balance the weight of the stone Mg in order for it to 
bounce. 



B. 



A circular stone 



For a circular stone, the immersed area is a more com- 
plex function of the height z, and is given in terms of the 
area of a truncated circle. A simple integral calculation 
yields 



S im (s) = i? 2 [arccos(l-s/i?)-(l-s/^)\/l - (1 - s/R) 2 } , 

(9) 

with s = | z | / sin 6 (the maximum immersed length) and 
R = a/2 as the radius of the stone. 



The equation of motion for z, Eq. (2b), thus becomes 
non-linear. However, it is possible to describe (at least 
qualitatively) the collisional process and obtain the con- 
dition for the stone to bounce. 

I first introduce dimensionlcss variables to simplify 
the calculations. The dimcnsionless height, z, time, r, 
and immersed area, A, are defined as z — —z/Rsm9, 
t = u)ot, and A(z) — S- lm /R 2 . (The minus sign in z is 
introduced for convenience.) If we use these variables, 
Eq. (2b), and V z = dz/dt, we obtain 



d 2 z 1 

a- -A{z) 



(10) 



with a = g/(Ru>Q sin 6). Equation ( [To| ) is the equation 
of a particle (with unit mass) in the potential V(z) = 
J (^A(z) — a)dz. We can use standard techniques for me- 
chanical systems to solve Eq. (^o|). In particular, Eq. (|Io| ) 
can be integrated once to give the "constant energy" con- 
dition 



+ V(z) = E, 



(11) 



where E is the energy of the system and is given in terms 
of the initial conditions 



1 



dz \ 2 



E =2\Tr 



T = 



V(5 = 0) = ~ (V z0 /(Ruo sin 8)) 2 . 

(12) 



The potential V(z) can be calculated analytically using 
the expression for the immersed area Si m given in Eq. (|^). 
An integral calculation gives 



2 1, s? 
- + -(1 - z) 2 

3 3 V ' 



-(1 — z)arccos(l — z) ] — az . 



(13) 



This potential is plotted in Fig. ||as a function of z. As a 
consequence of the constant energy condition, Eq. ([Tl]), 
z exhibits a turning point at a maximum depth defined 
by V(z max ) = E. 



V(z) 




FIG. 2: Plot of the potential V(z). The horizontal line is the 
constant energy E of the system. 

Here again, the condition for the stone to bounce is 
that this maximum depth be reached before the stone is 
fully immersed, that is, |z max | < 2i?sin#. In terms of di- 
mensionless variables, we obtain the condition: z max < 2, 
with z max defined by V(z max ) = \{V z0 /(Rujq sin6»)) 2 . 
This condition can be explicitly solved. Let me intro- 
duce zq such that dV /dz = at z — zq: Viz) is a mono- 
tonically increasing function of z for z > zq. Now it 
is easy to show that z max > zq (because V(z max ) > 
and V(zq) < 0), and the condition z max < 2 is there- 
fore equivalent to V(z max ) < V(2) = | — 2a, that is, 
±{V z0 /{Ru a sm6)) 2 < | - 2g/(Rw 2 sm9). Then the 
condition for skimming can be rewritten (recalling that 
V z0 /V x0 = tan/3) 



V x0 > V c 



16Mg 
7rCp w a 2 



l - 



8M tan 2 g 
7ra 3 Cp w sin t 



(14) 



Up to (slightly different) numerical factors this condition 
is the same as in Eq. (0) for a square stone. Note more- 
over, that the reasoning used for the potential V is quite 
general and can be applied to the square shape as well. 
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This reasoning yields the same condition as Eq. (g) in 
this case. 

Note also that for the circular stone, a simplified analy- 
sis of the motion could have been performed. First if z re- 
mains small during the bounce of the stone, a small z ex- 
pansion of V(z) is possible, yielding V(z) = 4y/2/15z 5 / 2 — 
az (corresponding to a parabolic approximation for the 
shape of the stone near its edge). Moreover, we remark 
that for small V z o, the energy E goes to zero, so that z max 
is defined in this case by V(z max ) = 0. If also use the pre- 
vious approximation, we obtain z max = (15a/4v / 2) 2 ^ 3 - 
The condition for the stone to bounce, z max < 2, there- 
fore yields a < 16/15. In terms of V x q, this condition 
gives again a minimum critical velocity for skimming, 
defined as V c = ^(Mg/Cp^ 2 with < = 15/4 ~ 3.75. 
This result is thus close to the "exact" condition found 
in Eq. (Oh for the V z0 = case. 



C. Energy Dissipation 

I have so far described the rebound of the stone by ana- 
lyzing its vertical motion. This analysis gave a minimum 
velocity for skimming which results from the balance be- 
tween the weight of the stone and the lift of the force due 
to water. However, some energy is dissipated during the 
collision due to the "friction" contribution of the force 
(the component along x) . This mechanism of dissipation 
leads to another minimum velocity condition, in terms 
of the balance between dissipation and initial kinetic en- 
ergy. Only a qualitative description of the dissipation is 
given here. 

As shown by Eq. (|2|), the component F x of the reac- 
tion force in the x direction (parallel to the water sur- 
face) will decrease the velocity of the stone. Then after 
a few bounces, the condition for the stone to bounce, 
Eq. (g) or Eq. (^i[), will no longer be satisfied and the 
stone will stop. It is possible to estimate the decrease 
in the x-component of the velocity using the equation 
for the center of mass position, Eq. (paj). If we multiply 
both sides of Eq. (|2a|) by V x and integrate over a collision 
time, we obtain the decrease in the kinetic energy in the 
x direction in terms of the work of the reaction force: 



1 



:MV„ 



1 



xf 



:MV X % = 



F x (t)V x (t)dt , (15) 



where V x q and V x f are the ^-components of the velocity 
before and after the collision, t co \\ is the collision time, 
and F x = \C p w V 2 S- lnl is the x-component of the reaction 
force, with C = Ci sin 9 + Cf cos 9. 

A rough estimate of the right-hand side of Eq. (|l5|) is 



F x (t)V x (t)dt~V x0 



toll 



F x (t)dt. (16) 



Now we have the simple relation F x (t) = fiF z (t), with 
/i = C/C (see Eq. (P). Moreover, it is expected that 
the average vertical force during a collision, (F z (t)) — 



t~ jj /o co " Fz(t)dt, is the order of the weight of the stone, 
M g. This point can be explicitly verified for the square 
stone case, using the expression of the force F z in terms 
of the height zlt) and Eq. (^). The final result is 
(F x (t)) ~ pA-Ig.u. Moreover, as shown in the above (and 
in particular for the square stone, although the results 
remain qualitatively valid for the circular one), the colli- 
sion time is given approximatively by t co u ~ 2tt /luq. We 
eventually find that the loss in kinetic energy in Eq. ( |l5| ) 
is approximatively given by 



W ~ -nMgV x0 — = -nMg£, 

UJq 



where £ is defined as 



„ 2tt „ /2A/sm<9 

V Cp w a 



(17) 



(18) 



The quantity i = V x t co \\ is the distance along x traversed 
by the stone during a collision. If the energy loss W is 
larger than the initial kinetic energy, the stone would be 
stopped during the collision. Using Eq. (|l5|), this condi- 
tion can be written explicitly as \MV x0 > \W\ = pMgi. 
We deduce that the initial velocity should be larger than 
the minimum velocity V c in order to perform at least one 
bounce, that is. 



v x0 >v c = ■ 



(19) 



If we use the same numerical values as in the previ- 
ous paragraph, we obtain p, — 1.4, £ = 13 cm, so that 
V c ~ 2ms~ 1 . This criterion is more restrictive than the 
previous one, Eq. (0). I thus consider in the following 
that Eq. (19) is the criterion for the stone to skim over 
water. 



V. WHY GIVE THE STONE A SPIN? 

The previous calculations assumed a constant angle 9. 
It is obvious that the rebound of the stone is optimized 
when 9 is small and positive (see, for example, the value 
of the force constant C = Ci cos 9 — Cf sin 9 which de- 
creases when 9 increases). Now, if after a collision, the 
stone is put in rotation around the y-axis (see Fig. [j]), 
that is, 9 ^= 0, its orientation would change by an appre- 
ciable amount during free flight: the incidence angle 6 
for the next collision has little chance to still be in a fa- 
vorable situation. The stone performs, say, at most one 
or two more collisions. There is therefore a need for a 
stabilizing angular motion. This is the role of the spin of 
the stone. 

Let us denote </>o as the rotational velocity of the stone 
around the symmetry axis parallel to n in Fig. |l|. I ne- 
glect in the following any frictional torque on the stone 
(associated with rotational motion). During the collision, 
the reaction force due to the water is applied only to the 
immersed part of the stone and results in a torque applied 
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on the stone. For simplicity, I consider only the lift part 
of the force. Its contribution to the torque (calculated at 
the center O of the stone) can be readily calculated as 
■A^iift = OP ■ FiiftBy, where e y is the unit vector in the y 
direction in Fig. [l] and P, the point of application of the 
lift force, is located at the center of mass of the immersed 
area. This torque is in the y direction and will eventually 
affect the angular motion along 9. However a spin motion 
around n induces a stabilizing torque: this is the well- 
known gyroscopic effect!] The derivation of the equation 
of motion of the rotating object (the Euler equations) 
is a classic problem and is treated in standard mechan- 
ics textbooks (see for example, Ref. [?]). On the basis of 
these equations, it is possible to derive the stabilizing gy- 
roscopic effect. This derivation is briefly summarized in 
the Appendix. 

In our case, the equation for the angle 9 can be written 



9 + lu 2 {9-9 ) 



Me 
Ji ' 



(20) 



where u> = [(Jo — Ji)/Ji]4>o, 4>q is the initial spin angu- 
lar velocity (in the n direction), and Jo and J\ are mo- 
ments of inertia in the n and t directions, respectively; 
9q is the initial tilt angle and Me = OP ■ Fua is the 
projection of the torque due to the water flow in the y 
direction. Equation (|2(i| ) shows that in the absence of 
spin motion, (f>o = 0, the torque due to the lift force will 
initiate rotational motion of the stone in the 9 direction. 
As discussed above, the corresponding situation is unsta- 
ble. On the other hand, spin motion induces a stabilizing 
torque that can maintain 9 around its initial value. The 
effect of the torque can be neglected if, after a collision 
with the water, the maximum amplitude of the motion 
of the angle 9 is small: 59 = [9 — 6* ] m ax «1. If we use 
Eq. (^0|) ) an estimate of 89 can be obtained by balancing 
the last two terms in Eq. (pfj|), yielding 89 ~ Me / '(J\uj 2 ) 
(note that up to numerical factors (Jo — Ji)/Ji) ~ 1 and 
J i ~ MR 2 , with R the radius of the stone). The order 
of m agnit ude of Me can be obtained using the results of 
Sec. IV C . The average vertical force acting on the stone 
has been found to be the order of the weight of the stone 
(see the discussion after Eq. Jl6[)): (F z (t)) ~ Mg. If we 
take OP ~ R, we obtain the simple result Me ~ MgR. 
The estimate for 89 follows directly as 89 ~ g/(Rio 2 ). 
Therefore, the condition for 9 to remain approximately 
constant, 89 <C 1, is 



to > 



(21) 



For a stone with a diameter of 10 cm, Eq. (|2l| ) gives 
<f)o ^> 14 s -1 , corresponding to a rotational frequency 
larger than a few revolutions per second (~ 2 Hz). This 
condition is easily fullfilled in practice and corresponds 
approximately to what we would expect intuitively for a 
successful throw. Note that the condition (|l]) is inde- 
pendent of the center of mass velocity of the stone V. 



VI. AN ESTIMATE FOR THE MAXIMUM 
NUMBER OF BOUNCES 

The estimation of the maximum number of bounces is 
the most difficult and tentative part of the analysis be- 
cause many factors can in principle slow down or desta- 
bilize the stone, some of which arc extremely difficult to 
model (such as irregularities of the water surface and the 
wind). We shall assume the idealized situation described 
above (perfect surface, no wind, idealized reaction force) 
and focus on two specific factors, which appear, at least 
intuitively, as natural candidates for stopping the stone. 



Slow down of the stone 



As I have discussed in Sec. [VC, energy is dissipated 
during a collision and the x component of the velocity of 
the stone will decrease during each collision: after a few 
collisions, all the initial kinetic energy will be dissipated. 
This process can be easily formulated. 

I consider a succession of N collisions. Between two 
collisions, the motion is parabolic (wind and air friction 
are neglected) and the initial x component of the velocity 
at the next collision is equal to the final x component of 
the velocity at the end of the previous collision. The 
important point to note is that the energy loss during 
one collision, Eq. (|l7|), is independent of the velocity V x o 
before the collision. Therefore, the velocity of the stone 
after N collisions obeys the relation 



1 



1 



MV 2 [N] - -MV 2 [0] = -NfiMgi, 



(22) 



so that the stone will be stopped at a collision number 
N c such that the total energy loss is larger than the ini- 
tial kinetic energy (similar to the argument leading to 
the critical velocity for skimming, V c , in Eq. (p"9| )). This 
criterion corresponds to V 2 [N c ] = in Eq. Q22|), and N c 
is given accordingly by 



Nr. 



(23) 



If we use the same typical values as before (M = 0.1kg, 
a = 0.1m, Ci « C f « 1, p w = lOOOnT 3 , /3 ~ ~ 10°, 
we obtain /i ~ 1.4 and i ~ 13 cm. We then find N c 6 
for the initial velocity V x q = 5ms" 1 , N c ~ 17 for V x q — 



ms , and N c w 38 for V x o = 12 ms 



The latter 



number of bounces corresponds to the world recordE 

It is interesting to calculate the distance between two 
successive collisions. As noted, the motion of the stone 
is parabolic out of the water: the position {X, Z} of 
the particle is given by X(t) — V x t, Z(t) = —\gt 2 + 
\V z \t. The next collision will occur at a distance AX = 
2V x \V z \/g. The dependence of V x on the number of colli- 
sions N is given by Eq. (p2|). On the other hand, V z does 
not depend on the number of collisions because the stone 
rebounds "elastically" in the z direction, as follows from 
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the analysis of the collisional process in Sec. [V (see, for 9 



example, the conservation of the energy E during the col 
lision discussed for the circular stone). If we use Eq. (22) 
we obtain the simple result 



AX[N] = AX Jl 



N 



(24) 



where AXo — 2V x o\V z o\/g. Note that AXq is approxi- 
mately equal to the distance between the two first rico- 
chets, AX[N = 1], when iV c > 1. For V x0 = 8ms" 1 , we 
obtain AX w 2.25 m. 

Equation @ for AX[N] is plotted in Fig. |. We re- 
mark that the decrease in the distance between two suc- 
cessive ricochets is first rather slow (AX[N) ~ AX (1 — 
2^-) for N <C N c , see Eq. (|I|)), but strongly accelerates 
for the last collisions when N ~ N c , due to the square 
root variation of AX[N] close to N c . This result is in 
agreement with observation. Such an effect, is known to 
specialists of stone-skipping as "pitty-pat."El 



is 69 ~ g/(Ru> ), with lu ~ (f>o, the (constant) spin 
velocity of the stone. Now assume that the destabiliz- 
ing effects add, a reasonable assumption. Then, after N 
collisions we expect that A^9 ~ N59. The stone is com- 
pletely destabilized for a collision number N c such that 
A N J ~ 1, yielding 
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(25) 



If we use the same numerical values as before, we ob- 
tain, for example, N c ~ 5 for a initial spin velocity 
4>q = 5rev/s and N c — 38 (the world recorco) for 
4>o = 14rev/s. Note, however, that there is a quite large 
uncertainty of the numerical prefactors in the above es- 
timate of N c , and this estimate is merely qualitative and 
should not be taken literally. 



VII. DISCUSSION 



x 



X 



0.5 




FIG. 3: Plot of the (normalized) distance between two suc- 
cessive collisions AX [N]/ AXo as a function of the number of 
bounces N. The initial velocity is V x o = 8 ms -1 , correspond- 
ing to N c = 17 (using the same values for the parameters as 
those given in the text). The vertical dashed line indicates 
that N c = 17. 



B. Angular destabilization 

However, there is another possible destabilizing mech- 
anism in the collision process. As was discussed in Sec. 0, 
the rotational stability of the stone is crucial in the colli- 
sional process. A criterion for stability has been found in 
the form of a minimum spin velocity of the stone. How- 
ever, each collision will perturb the rotational motion and 
the sum of all these effects can eventually bypass the sta- 
bility condition. This argument can be easily formulated. 
As shown above, the amplitude of the angular motion of 



At the level of our description, the maximum number 
of bounces results from the combination of the two previ- 
ous mechanisms: slow down and angular destabilization. 
The maximum number of bounces is therefore give n by 
the minimum of the two previous estimates, in Eqs. (|23|) 
and @. 

The estimate N^ d obtained in Eq. (^j|) from the slow 
down of the stone depends only (quadratically) on the ini- 
tial velocity of the stone: in principle, a very large num- 
ber of bounces could be reached by increasing the initial 
velocity of the stone. But on the other hand, the angu- 
lar destabilization process results in a maximum value of 
jyspin w i]j c h is independent of the initial velocity of the 
stone, as indicated by Eq. (25). This shows that even 



if the initial velocity of the stone is very large, that is, 
A^ d ^> 1, the stone will be stopped by angular destabi- 
lization after 7V^ pin bounces. In other words, the initial 
"kick" that puts the stone in rotational motion is a key 
factor for a good throw. 

The results presented here are in agreement with our 
intuition for the conditions of a good throw. Some of the 
results are also in agreement with observations, for ex- 
ample, the acceleration of the number of collisions at the 
end of the throw (a phenomenon known as "pitty-pat" in 
stone skipping competitions!!!) . Some easy checks of the 
assumptions underlying our calculations could be per- 
formed, even without any sophisticated apparatus. For 
example, taking pictures of the water surface after the 
ricochets would locate the positions of the collisions (be- 
cause small waves are produced at the surface of water). 
A simple test of the variation of the distance between 
two collisions as a function of collision number, Eq. (B3), 
would then be possible. A more ambitious project would 
be to design a "catapult," allowing one to throw stones 
with a controlled translational and spin velocity (together 
with the incidence angle of the stone on water) . A mea- 
surement of the maximum number of bounces performed 
for various throw parameters would allow us to check the 
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assumptions underlying the present simple analysis and 
to determine some of the parameters involved in the de- 
scription (such as \x and t) . It would be also interesting to 
repeat the experiments reported in Ref. |^ using modern 
techniques (such as fast cameras), in order to image and 
analyze in particular the rebound process as a function of 
the throw parameters. Hopefully a better understanding 
of the mechanisms of stone skipping will allow someone 
to break the actual world record! 
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APPENDIX A: APPENDIX 

I briefly recall the derivation of Eq. ( |2(i| ) , from the Euler 
equations described in Ref. The latter are written as:EI 



j duj\ 

J l rlt 



dt 
t duj 2 



■ u 2 ui 3 (h - I 3 ) =Ni (Ala) 
uHJaih-h) = N 2 (Alb) 
uo^ih-h) =N 3 . (Ale) 



In Eq. (Al), I a , uj a , and N a (a — 1, 2, 3) are respectively 
the moment of inertia, angular velocity, and torque along 
the direction of a particular principal axis, denoted as 
a. In our case, the direction 1 is taken along the axis 
perpendicular to the vectors n and t (the direction 1 is 
along the y axis in Fig. [l]) , the direction 2 along n and the 
direction 3 along t. We therefore have ui± = 9, and due 
to the symmetry of the circular stone, I\ = 1% = J\ and 
I 2 = Jo- Moreover, because only the lift component of 
the reaction force (along n) is considered in the present 
analysis, we have Ni = M.$ and N 2 = N% = 0. 

Equation (All]) yields immediately that ui 2 = 0. We 
therefore have w 2 — cf>o, with cf>o the initial spin velocity. 



Equation (Ale) can be therefore written as: 



duj3 
dt 



J i — ,Jq 
Ji 



(A2) 



If we use tui — 9, Eq. (|aJ) can be integrated once to give 



w 3 



Ji — Jo 



(A3) 



with 9p 9 (t = 0), th e init ial tilt angle. The substitution 

of Eq. (|A|) into Eq. (|Ala|) leads to Eq. @). 



The actual world record appears to be 38 re- 
bounds (by J. Coleman-McGhe e) . See, for example, 
<tittp : //www. stoneskipping. com> for more information 



on stone skipping competitions. 
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